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I. INTRODUCTION 

111 the first part of our paper [ij, extending an approacli formulated in Ref. to tlie relativistic limit, we have 
studied systematically the structure of the singularities formed in a relativistically large amplitude plasma wave close 
to the wavebreaking in a thermal plasma. We have shown that typically the electron density distribution in the 
breaking wave has a "peakon" form with a discontinuous coordinate dependence of its first derivative, similar to the 
profiles of nonlinear water waves and that in the above breaking limit the derivative becomes infinite. This 

results in a finite reflectivity of an electromagnetic wave interacting with nonlinear plasma waves. In particular, this 
is an important property because nonlinear Langmuir waves play a key role in the "relativistic flying mirror" concept 
(^-[Hl- In this concept, very high density electron shells are formed in the nonlinear wake wave generated by an 
ultrashort laser pulse propagating in an underdense plasma with a speed close to the speed of light in vacuum. The 
shells act as mirrors flying with relativistic velocity. When they reflect a counterpropagating electromagnetic pulse, the 
pulse is compressed, its frequency is upshifted and its intensity increased. It is the singularity in the electron density 
distribution that allows for a high efficiency in the reflection of a portion of the counterpropagating electromagnetic 
pulse. If the Langmuir wave is far below the wave-breaking threshold, its reflectivity is exponentially small. For a 
nonlinear Langmuir wave the singularity formed in the electron density breaks the geometric optics approximation 
and leads to a reflection coefficient that is not exponentially small @, ■ 

In the present paper we address the problem of the interaction of an electromagnetic wave with a nonlinear plasma 
wave which is of interest for the "photon accelerator" concept [T^ and for the "relativistic flying mirror" paradigm 
[6l-[lH . We calculate the reflection coefficients of an electromagnetic wave at the singularities of the electron density 
in the most typical regimes of a strongly nonlinear wave breaking in thermal plasmas. 



II. ELECTROMAGNETIC WAVE REFLECTION BY THE ELECTRON DENSITY MODULATED IN 

THE BREAKING WAVE 

As we have seen in the first part of our paper [l[ , in a strongly nonlinear wake wave the electron density is modulated 
and forms thin shells (singularities or caustics in the plasma flow) moving with velocity /3ph. In the Introduction, in 
a way of Refs. we have discussed how a counterpropagating electromagnetic wave can be partially reflected 

from these density shells which play the role of relativistic mirrors. While in the case of a cold plasma the electron 
density at the singularity tends to infinity (see Eq. (59) of Part I [l[ and Refs. 0, EH), in a thermal plasma the 
density is limited by the expressions given by Eqs. (41) and (42) of Part I Although in this case the density 
profile is described by a continuous function of the variable X, its derivatives with respect to X are discontinuous. 
This discontinuity results in the breaking of the geometric optics approximation and leads to a reflectivity that is not 
exponentially small. 

In order to calculate the reflection coefficient, we consider the interaction of an electromagnetic wave with the 
electron density shell formed at the breaking point of a Langmuir wave in a thermal plasma similarly to what has 
been done in Refs. The electromagnetic wave, described by the z component of the vector potential Az{x, y, t), 

evolves according to the linearized wave equation 

dttA, - {da,,A, + dyyA,) + ^l^{x - v^ht)A, = 0, (1) 

where we have reverted to dimensional units and 

The last term in the l.h.s. of Eq. ([5]) is the z— component of the electric current density generated by the elec- 
tromagnetic wave in a plasma with the electron distribution function fe{p)- In the limit 7phApo ^ 1 for the 
electromagnetic wave frequency larger than the Langmuir frequency calculated for the maximal electron density, 
Lo 3> a;pc(27ph/Apo)^''^, we can neglect the finite temperature effects on the electromagnetic wave dispersion, which 
have been analyzed in Ref. [l3l |. in the limit of homogeneous, stationary plasmas. For the water-bag distribution 
function 

/e(p, X) = n„e{p - p^{X))0{p+{X) - p)/Apo (3) 

ripg(X) takes the form 



p+ix) + ^i+p+{xy 

= Apo \p-{X) + Vl+P-(^)2 



Keix) = In ( i ^ (4) 
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where = innoe'^ /rUe, P±{X) and Apo are now dimensionless (normalized on niec). 

The wake wave modulates the electron density and temperature increasing them in the compression regions and 
decreasing them in the rarefaction regions. In Fig. [1] we illustrate the dependence of npe{X)/ujpe on X for the 
parameters of a wakewave corresponding to Apo = 0.1 and £',„ax = 2.3 at X = 15 and for /3ph = 0.992. 
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FIG. 1. Dependence of the frequency ratio Qpe/ujpe on the coordinate X for the parameters of a wakewave corresponding to 
Apo = 0.1 and iSmax ~ 2.3 at X = 15 and for /3ph — 0.992. In the inset the ratio ^lpe{X) /cupe is shown in the vicinity of the 
maximum. 

From Eqs. (35) and (59) of Part I [3| in the ultrarclativistic case, /3ph ~ 1, using Eq. (j4]) wc find for a relatively 
cold distribution such that p_ <C 1 that near the wavebreaking point Q,p^ {X) is given by 



^UX) . ^ - ^fc^lxl, (5) 

7ph Apo 

The propagation of a sufficiently short electromagnetic wave packet in the plasma with electron density modulated by 
the Langmuir wave can be described within the framework of the geometric optics approximation. The electromagnetic 
wave is represented as a particle ( "photon" ) with coordinate x and momentum k (wave vector) . The interaction of a 
"photon" with a Langmuir wave that propagates with a relativistic phase velocity Vph ~ c can be accompanied by a 
substantial frequency upshift called "photon acceleration" [l^ - [l7| . Using the dispersion equation 

w(x, k; t) = ^k^c^ + ni^{x~vpi,t), (6) 

where fc^ = fc^ + k'j_ with fc|| and k± the wave vector components parallel and perpendicular to the propagation 
direction of the Langmuir wave, we obtain the "photon" Hamiltonian function which depends on the canonical 
variables X = x — Vpht and fc|| (see Ref. [2^ ) 

Hphoton(^,fc||) = y'fc2c2+r!2g(X)~/3phfc||C. (7) 

The transverse component of the wave vector is constant k± = fc^ o ^md k±^Q = is assumed for the sake of simplicity. 

The phase portrait of the photon for il,pe{X) given by Eq. (|4| for the parameters corresponding to Fig. [l]is shown 
in Fig. m 

Along an orbit corresponding to the value 'Hphoton(-'^, ) ~ 'Hphoton(-''^o, fc||,o) ~ 'Hphoton^o of the Hamiltonian ([7]) 
the photon frequency is given by 



W — 7ph^photon,0 



l±/3phjl 



^photon, OTph 



(8) 



Photons, for which ?^photon,o < niax{r2pe/7ph} are trapped inside the region encircled by the separatrix. Along the 
orbit their frequency changes from Wmax and Wmin corresponding to the plus and minus signs in the r.h.s. of Eq. ([8]) 
at the minimum of Vlpe^X). 

Photons with 'Hphoton,o > ^8J^{^pe/lph} a-re not trapped and for them the sign in the r.h.s. of Eq. ([5]) does not 
change. For trajectories far above the separatrix the photon frequency variations are relatively weak. However a 
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FIG. 2. Photon phase portrait for the parameters of a nonhnear Langmuir wave corresponding to Fig. \T\ The dashed line 
corresponds to the trajectory of photons that have appeared due to the over-barrier reflection at the crest of the breaking wave. 
In the inset the photon trajectories in the vicinity of the saddle point are shown. 



sufficiently strong wakefield can reflect a counterpropagating photon, q < due to above-barrier reflection (this 
trajectory is shown in Fig. [5] by a dashed line). Such a photon acquires a frequency 

according to the Einstein formula for the frequency of the electromagnetic wave reflected by a relativistic mirror 
p^ . The geometric optics approximation fails when the wakefield is close to wave breaking and this provides the 
appropriate conditions for a not exponentially weak wave scattering. 



III. ABOVE-BARRIER SCATTERING OF AN ELECTROMAGNETIC WAVE AT THE CREST OF THE 

BREAKING WAKE WAVE 



In order to find the reflectivity of the nonlinear wake wave we make a Lorentz transformation to the frame of 
reference moving with the phase velocity of the Langmuir wave. In the boosted frame, Eq. ([S]) for the electromagnetic 
wave interacting with the nonlinear Langmuir wave can be written as 



with 



-^g^(C)a(C)=0 (10) 



eA 

a(C) = ^ exp [-i{uj't' - kyy)] (11) 

TOeC 

and in the neighbourhood of the breaking point (C) can be written as 

q\C) = s'+9-i\C\. (12) 

Here 

— —5 k'^ , (13) 

" c^7ph 

and C = ^7phj t', k' , uj' are the coordinate and time and the wave number and frequency in the boosted frame of 
reference. The coefficient g-i is equal to 

9-1 = . '"l (14) 



FIG. 3. Scattering geometry. 



We seek for the solution to the above-barrier scattering problem for Eq. pU)) writing its solution in the form (see 
Refs. [ilillll) 

a(C) = [6+ exp (iW(S,)) + h_ exp {-%W{C,))\ , (15) 



where the phase integral is defined as 



M^(C)=/%(C'K'- (16) 



The above-barrier scattering geometry is illustrated in Fig. [3l For constant 6+ and &_ Eq. HH) corresponds to 
the "WKB" solution |25|. In the following the coefficients 6+ and 6_ are considered as functions of W instead of C, 
because, as explained in Ref. [llj, the mapping between W and ^ given by Eq. ([T6| is one-to-one on the real axis. Far 
from the breaking point, i.e. formally for ^ — >■ ±oo the function q'^iC.) reduces to a constant and the solutions 

arc exact, so that b± — >■ const as W ~> ±oo. 

In other words, the boundary conditions at C — > ±oo are 

6+(+oo) = l, 6_(+oo)=p, , . 

5+(-oo) = 0, 6_(-cx))=r. 

Since in the representation (|15|) . the single unknown function a(Q has been replaced by the two unknown functions 
^±(C)i ^ subsidiary condition is necessary. We shall impose the condition 

|=^V^(M^'^(C)-6_e-^(^)). (18) 

Differentiating Eq. p3|) with respect to ( and taking into account the constraint p^ . we find 

dQ dQ dC, \ / 

while differentiating Eq. (|T8|) with respect to C and substituting (PajdC^ into Eq. (|10p yields 

^e'M/(C) _ ^^-^w(Q ^ / ^_^-w(C) _ i^^wiQ\ (20) 

dC (iC dC, \ / 



The system of Eqs. ((T^ and (PH)) is equivalent to Eq. pUj) . It can be rewritten in the form 

_d_( h+ \ _( 5(VF)e2»^ W 6+ 

dW\b- ~\ 5(W^)e-2'^ j I 6_ 



(21) 
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with 



S{W) = l^lnq{aW)). (22) 



For q{() given by Eq. ([T2|) we have 



WiO = —- \{s'+g-iK\f' - s'] sign(C), (23) 



where sign(C) = — 1 if C < and sign(C) = 1 for ^ > 0, and 



/o \ 1/3 

9(1^) = 4r^M^sign(C) + s3 (24) 



so that 



^^^^ = i?f(W^^^"^^^- ^''^ 
It follows that S{W) is discontinuous at C (W 0) 

5(VK = 0) = gsign(C). (26) 

Integrating both sides of Eq. (1211) and using the above formulated boundary conditions for 6-|-(±oo), we can obtain 
the reflection coefficient p in the form of the infinite series (26j 

oo nWn-1 r + oo 

P = -5](-l)™/ dW^o5(W^o)e2*'^''n / dVnS{Vn)e-^'''- dW^„5(W^„)e2'^", (27) 

with the product equal to unity for m = 0. 

The function q{() defined by Eq. (|12p has a discountinous first derivative at C = 0. In the vicinity of the singularity 
point it can be represented in the form q{Q w go + ^ZilCI with qo ~ s and qi = g_i/2s. Expanding W{C) and S{W) 
in powers of ^ and substituting them into Eq. ([27l) we can find (see Eq. (27) of Ref. [2l|) that the first term yields 
the dominant contribution to the reflection coefficient, with the result 

-tqi _ -^g-l 

and 

i?_i = |p_iP=5-i^- (29) 

Applicability of the WKB theory implies that f/_i ^ s'^. 

Similarly (see also Ref. [l2j) we can find the refiection coefficient at the electron density singularity formed in the 
above breaking regime discussed in Part I In this case the electron density distribution is given by Eq. (106) of 
Part I. Using this relationship we obtain 



/ + 00 
exp [2z<] (^(OyC - ^(C - AOTC^j dC (30) 



where 



kp = c/ujpe and = Apo/ei5max- Calculating the integral pop we find 



,3/2 3/2 V^efimaxmeC 

%4)=V 7ph , (31) 



n I ■^ /T^ exp (zsAQ sin (sAQ , , 
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Consequently, we write 



R 



An ■ 



(33) 



From Eqs. ((29)) and ([33]) we can see that in thermal plasmas the reflection coefficient is s ^ 1 times larger in the 
above breaking regime than for a wake wave approaching the wavebreaking threshold. 

Generalizing Eqs. (106) and (108) of Part I [l[, we can write the electron density dependence on the coordinate ( 
in the form 



«e(C)-^[e(C+)(C+)^/™-^?(C-)(C-)^/^ 



(34) 



with C± = C T and m an even number, and 



neiC) 



no 
AC 



' (C+) (C+)'/" + (-C-) (-C-)'/™ -oiC-) iC-f"- - e (-C+) (-C+)'/" 



(35) 



for m an odd number. 

It is easy to show that for the reflection coefficient, j^-j = '^6 have 



4(-^.)-v™r ( 1 + 1^ "^^(-^^^^ 



(36) 



if m is even, and 



r)=^'(^,^) 



2(i + (-i)i/'")(-is)-i/™r ( 1 + 



1 \ sin(sAC) 



sAC 



(37) 



if m is odd. 

Comparing Eqs (|331 [36l [37)) for the reflection coefficient with the corresponding coefficients obtained in Ref. [l^l , 
we find that the effects of a finate temperature enter Eqs. ([33l [36l [37| as a form-factor |sin(sAC)/sACp. In the limit 
ApQ this form factor tends to unity while for sAC, ^ 1 decreases as « l/(sAC)^. 

Since the frequency, uJr, and the number of refiected photons, N^, are related to that incident on the relativistic 
mirror wq and Ns as ujr = cjo(l + /3ph)/(l — Pph) ~ '^o47ph and Nr = RNs, the energy of the refected photon beam 



is given by £r ~ ^s^lph-^^ where £5 is the energy of the laser pulse incident on the mirror. Comparing with the 
energy of the electrons in the first period of the wake wave (e.g. see [111), ~ '^^ias,d('^pc/wo)^, where fias,d is the laser 
driver energy, we find that the photon back reaction (the ponderomotive pressure) on the wake wave can be neglected 
provided £3 ^ '^ias.d/47piji?. As a typical reflection coefficient value we can take R ~ l/7ph (s^e Refs.[i,[i|) and 
obtain the condition of relative weakness of the incident laser pulse £3 < fias,d- As we see owing to the weakness of 
the photon-wakewave interaction the incident pulse energy can be of the order of that in the driver laser pulse. 



IV. DISCUSSIONS AND CONCLUSIONS 



In the first Part of our paper [l| we found the structure of the typical singularities that appear in the electron density 
during the wave breaking in a thermal plasma. The singularity in the electron density, moving along with the wake 
wave excited by a high intensity ultra-short pulse laser, can act as a fiying relativistic mirror for counterpropagating 
electromagnetic radiation, leading to coherent reflection accompanied by the upshift of the radiation frequency. This 
process implies finite (not exponentially small) refiectivity at the electron density singularities. This is provided by the 
structure of the singularity formed in a relativistically large amplitude plasma wave close to the wavebreaking limit 
that leads to a refraction coefficient with discontinuous coordinate derivatives. We found the reflection coefficients of 
an electromagnetic wave at the singularities of the electron density in the most typical regimes of strongly nonlinear 
wave breaking in thermal plasmas. The efficiency of the photon refiection can be substantially increased by using the 
above breaking limit regimes which lead to the formation of high-order singularities. 
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